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Periodic structures with subwavelength features are instrumental in the versatile and effective control of elec- 
tromagnetic waves from radio frequencies up to optics. In this paper, we theoretically evaluate the potential ap- 
plications and performance of electromagnetic metasurfaces made of periodically patterned graphene. Several 
graphene metasurfaces are presented, thereby demonstrating that such ultrathin surfaces can be used to dynam- 
ically control the electromagnetic wave reflection, absorption, or polarization. Indeed, owing to the graphene 
properties, the structure performance in terms of resonance frequencies and bandwidths changes with the varia- 
tion of electrostatic bias fields. To demonstrate the applicability of the concept at different frequency ranges, the 
examples provided range from microwave to infrared, corresponding to graphene features with length-scales of 
a few millimeters down to about a micrometer, respectively. The results are obtained using a full-vector semi- 
analytical numerical technique developed to accurately model the graphene-based multilayer periodic structures 
under study. 

PACS numbers: 33.57,+c and 81.05.Xj 



I. INTRODUCTION 

Graphene is a flat monolayer of carbon atoms arranged in a 
honeycomb lattice, which can also be viewed as the building 
block for other graphitic materials, such as fullerenes, carbon 
nanotubes and bulk graphite^. The experiments on graphene 
confirmed that its electrons behave as massless Dirac fermions 
with the possibility of traversing long distances without being 
scattered by the crystal, which in turn results in the capabil- 
ity to sustain large electric currents 3 6 . These unprecedented 
properties underpinned extensive applications and develop- 
ment of ultrathin carbon nanoelectronic active devices using 
small graphene samples^. In the early years of graphene re- 
search, it was one of the most expensive materials obtained 
through the nanofabrication techniques 9 . However, recently 
efficient growth of large graphene samples has been achieved 
using chemical vapor deposition, which has now become a 
common graphene fabrication techniqu e 1 ™ 11 ! This allowed 
the production of larger graphene based structures as needed 
to realize the graphene electromagnetic (EM) metasurfaces 
discussed in this paper. 

One novel area, where the properties of graphene may 
be influential, is the dynamic control of the EM waves 
propagatiorfGK]. Notably, the capability of graphene to serve 
as a platform in transformation optics has been shown based 
on a homogeneous (i.e. unpatterned) but spatially bias mod- 
ulated graphene layer represent ed by a free-standing sur- 
face in a commercial EM soiverEIEl More recently, em- 
ploying graphene to realize cloaking surfaces has also been 
proposed-^. Furthermore, employing a graphene monolayer 
for sensing and measuring a metamaterial performance was 
also suggested 13 . These studies on graphene EM metasurfaces 
considered a homogeneous layer of graphene. Only very re- 
cently, research effort has been devoted on nanostructuring 
graphene to control its effective properties. For instance, the 
use of graphene strips for the realization of plasmonic waveg- 
uides with confined field profiles have been investigated 17 . 
The possibility of controlling the plasmonic resonances in a 



periodic arrangement of such strips has also been reported 18 . 
The effect of patterning on the transmission properties of this 
monoatomic layers were studied as welP. Finally, the use of 
periodic graphene metasurfaces to achieve strong absorption 
of infra-red radiation was reported 2 ^" 2 ^. 

The biperiodic patterning of metal has been long used at mi- 
cro and millimeter-wave frequencies to achieve electromag- 
netic properties which cannot be obtained using uniform mul- 
tilayered structures. Perhaps, the most widespread application 
of such periodic screens lies in frequency-selective surfaces 
(FSS), which are spectral filters for free space propagation 
applications (as employed in radar technologyj 2 ^ 2 ^! Though 
different FSSs having fixed frequency responses are success- 
fully modeled, implemented and utilized, their dynamic con- 
trol is still a challenging issue. The mechanism allowing 
the dynamic control of the FSS performance generally in- 
troduces very tough limitations on the achievable character- 
istics. In the microwave regime, the use of liquid crystals, 
varactor diodes, and MEMS switches t o p rovide a dynamic 
frequency response has been reportecPSEir However, achiev- 
ing such characteristics in the terahertz and infrared regime is 
even more difficult because of the aforementioned technolog- 
ical reasons. In this context, the use of micro/nano-patterned 
graphene seems extremely interesting to overcome the limita- 
tions of existing technologies, in terms of operation frequency, 
biasing c omple xity (thanks to graphene well-known electric 
field effec t 30 ! 3 ' I), as well as integration and miniaturization. 

We demonstrate this potential by presenting the model- 
ing and design of graphene-based dynamically-controllable 
metasurface, at both microwave and THz frequencies. First, 
graphene conductivity model is recalled and the numerical al- 
gorithm used for modeling the electromagnetic interaction of 
a plane wave with a patterned graphene metasurface is pre- 
sented. Subsequently, application examples are analyzed us- 
ing the proposed method. It is shown that electronic gating 
of a miniaturized graphene layer allows one to control both 
the appearance and the position of the periodic structure reso- 
nances. Results for a controllable X-band microwave absorber 
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and an infrared switchable polarizer are presented, highlight- 
ing the possibility to also affect absorption and polarization 
through anisotropic structuring. The designed structures are 
based on multilayer graphene stacks, since conventional bi- 
asing techniques cannot be employed without shielding the 
metasurface effective response (note that a similar idea has re- 
cently been proposed for fabricating THz modulators^). The 
different examples show that graphene not only enables con- 
venient dynamic control of the structure response, but also that 
its periodic patterning also allows tailoring the response both 
in terms of general functionality and operation frequencies. 



has been the subject of many recent research efforts, includ- 
ing successful comparisons with experimental resultP^JEH 
However, like any other material properties, these coefficients 
are strongly influenced by the fabrication process and envi- 
ronmental effects. For instance, the defects introduced by the 
substrate in the graphene sheet often incur deviations from the 
theoretical predictions 39 . In the presented study, the numeri- 
cal routine is developed for general parameters and later in 
the numerical examples values obtained based on theoretical 
models are considered. 

The conductivity value a can be theoretically calculated us- 
ing Kubo's formalism 30 , which yields the following equation: 



n. GRAPHENE CONDUCTIVITY MODEL 

A graphene sheet can be modeled as an infinitesimally 
thin conductive sheet, whose conductivity is obtained using a 
semi-classical quantum mechanical method. In low frequency 
regime (low THz regime and below) and in the absence of 
magnetostatic biasing, this conductivity is represented by a 
scalar value 33 . However, at higher THz and mid infra-red fre- 
quencies for very intense spatial variations of the interacting 
wave (in the scale of the scattering length), the conductiv- 
ity may possess non-negligible non-diagonal terms due to the 
spatial dispersion effec t 33 ^ 34 !. Additionally, in the presence of 
a magnetic bias field, a similar dyadic conductivity is needed 
to model the graphene sheet. Therefore, the conductivity is 
generally modeled as the following tensor: 
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where M n = \J A 2 + 2nvp\eBo\h, and f d (e) is the Fermi - 
Dirac distribution given by 
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where io is the radial frequency, [i c is the chemical potential 
hinging upon the applied electrostatic bias field Eo = zEq 
or doping, T is a phenomenological electron scattering rate, 
T is the temperature and Bo = zBq is the applied magneto- 
static bias field. In this study, we consider graphene periodic 
surfaces under electric bias field only, thus the anisotropic ef- 
fect emanating from magnetic biasing is set to zero. The four 
elements of the conductivity tensor can be written in the fol- 
lowing general fonrl 3 - 3 ^ 

d 2 n d 2 

dy 2 



'xy 



'yy 



2/3 
2/3 



dx 2 
d 2 



dxdy 

d 2 



dxdy 

d 2 d 2 



dx 2 



dy 2 



(2) 



The operator terms in |2]) represent the spatial dispersion ef- 
fect and are negligible either in low frequency regime or high 
phase velocities. 

The coefficients in namely a, j3, and a, can be obtained 
experimentally or analytically based on the existing models. 
In fact, the theoretical modeling of graphene conductivity 
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The above equation returns the conductivity for a general case, 
with both electric and magnetic bias fields. However, for the 
case of no magnetic bias field, the above equation should be 
calculated when |e£?o| tends to zero. In this case, the follow- 
ing equation should be usecPJ 
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Room temperature (T = 300 K) is considered throughout the 
paper. Correspondingly, the excitonic energy gap A is set to 
zero, the Fermi velocity in graphene is vf = 10 6 m/s and V is 
set to 12.2 meV 6 l 39 l 41 l Note that for the case of no magnetic 
bias field, the above equation should be calculated when \ eBo\ 
tends to zero. Fig. [T] shows the variation of graphene conduc- 
tivity with frequency for different bias electrostatic fields in 
the microwave and infrared regime. In the microwave regime, 
though the conductivity does not change significantly with 
frequency, it is strongly dependent on the bias electric field. 
The curves evidence stronger sensitivity of the conductivity to 
the bias field at low THz frequencies compared to the higher 
ones. It can be observed that changing the bias electrostatic 
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FIG. 1: (a) Real part and (b) imaginary part of the conductivity (ad) 
in terms of frequency for various bias electric field in microwave 
regime 



field values directly affects the range of frequency over which 
graphene behaves as a conductor. 

The coefficients a and f3 are obtained using a perturbation 
theory approach presented in 33 . Although this was initially 
made for unbiased graphene, following the same approach it 
can be shown that the final results also hold for electrically 
biased graphene. The equations providing the values of a and 
/3 then reads as 
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III. PMOM FOR GRAPHENE BIPERIODIC SURFACES 

Modeling EM properties of homogeneous grap hene lay- 
ers has been targeted in some previous publication d 33 ! 40 ! 42 * 4 ^, 
where the solution of the Maxwell equations in a one- 
dimensional space is used to simulate multilayer structures. 
The analysis of periodically patterned graphene demands 
more advanced analysis techniques, to account for the the 



periodic boundary condition on the electromagnetic fields as 
well as the boundary condition on the induced currents, i.e. 
normal component at the graphene boundaries should vanish. 
The properties of an EM screen can be obtained by formulat- 
ing the scattered field as a result of a plane wave excitation. 
The whole metasurface characteristics, such as resonance fre- 
quencies, high absorption frequency points, band diagram and 
polarization sensitivity, are then extracted from the scattered 
field spectrum. Here a periodic method of moments (PMoM) 
technique is implemented, which efficiently solves Maxwell's 
equations for multilayer planar metasurfaces 44 . This full- 
vector semi-analytical numerical technique not only allows 
accurately modeling periodic graphene surfaces but also in- 
cludes the effect of the substrate actual topology. 

The concept of PMoM for periodic surface s of va rious 
types has been introduced in many publication ;) 24 ! 25 ! 44 !. The 
method begins with the formulation of the boundary condi- 
tion which must hold on the graphene surface, i.e. 



K 

Ei 



Ef, 



(7) 



where Ef = [E s x ,E s y ] T and E* - [££ 



Ey) T are tangen- 



tial components of the scattered and incident electric field on 
graphene, respectively. In addition, Z s stands for the surface 
impedance and J = [J x , J y ] T is the induced current on the 
graphene surface. The scattered field can be written in terms 
of a Green's function and the induced currents. In case of a 
homogeneous substrate, this yields the following equation: 
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where G xx , G xy , G yx and G yy are the components of the 
dyadic Green's function in the spectral domain. k xm and k yn 
are given by 
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where ki nc = k x x + k y y is the wave vector of the incident 
plane wave and the pair (L x , L y ) determines the lattice con- 
stants of the periodic structure in both x and y directions, re- 
spectively. 

If the summations in ^ are cast in a matrix form, the fol- 
lowing equation is obtained: 
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(ii) 



[exp(j k xm x + jk yn y)] T is a row matrix containing the expo- 
nential terms and [0] T is a zero matrix with the same size as 
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and d/dy 



[exp(jk xm x + jk yn y)] T . In fact, the superscript T represent- 
ing the transpose sign is used to distinguish between row and 
column vectors. J x and 3 y are column vectors obtained from 
the Fourier coefficients J Xmn and J Vmn in (k xm ,k yn ) basis. spectral domain is as the following: 



-jk y should be used, where k x and k y are 



diagonal matrices with diagonal elements equal to k xm and 
k yn . Therefore, the conductivity equation for graphene in the 



Finally, G is the Green's function matrix. 



For solving (10 1 using the concept of the MoM, electric 
currents excited on the patches should be expanded in terms 
of some basis functions 
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where 3 X and J y are row vectors containing the basis func- 
tions used for expanding J x and J y , respectively. The un- 
known coefficients of these functions are arranged in the col- 
umn vector C. Using Galerkin's method and after some al- 
gebraic operations the following system of equations is ob- 
tained: 



[J X ] T [J,/ (G + Z 
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where J = 3 x x + J y y and E' = E x x + E y y is the incident 
electric field vector. [J x ] and [J y ] are matrices whose fc'th 
columns are Fourier coefficients of £'th corresponding basis 
functions. The signs * and t stand for the c omp lex and Her- 
mitian conjugate respectively. J x and J y in ( 10 1 are related to 
[J x ] and [J y ] through 
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Using the obtained coefficients C, all the desired quantities 
such as reflection and transmission coefficients can easily be 
calculated. The term exp(jki nc • r) is considered as a phase 
factor in all basis functions. Therefore, the Fourier coeffi- 



cients are calculated in ( 
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The studies on modeling planar geometries containing pe- 
riodic patches are mainly limited to considering various sub- 
strate properties. Different schemes are developed for homo- 
geneous, lossy, multilayer 2425 , periodic^ and anisotropic! 4 ^ 4 ^ 
substrates. Thin metallic patches are accurately modeled with 
an equivalent scalar surface impedance. Therefore, a scalar 
value is always assumed for the matrix Z s . However, for 
graphene patches, a more general form should be developed 
for Tig, which is actually the goal in this section. 

Let us first take the conductivity equation J t = erE t into 
account. By considering the operator form of the graphene 
conductivity, the general equation for graphene reads as 
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For the analysis of a frequency selective surface using PMoM, 
this equation should be transformed in the spectral domain. 
For this purpose, the equivalence equations d/dx = — jk x 
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By simply comparing ( [To} with (JTOj, the matrix Z s can be 
deduced as 



a — ak x — /3kj/ 
-2/3k x k y ' 



-2/?k x k y 
a — /3k x 2 — ak y 



(17) 



Note that the inverted matrix consists of four diagonal ma- 
trices. Thus, its inversion can be accomplished analytically, 
without the need to follow complicated computational proce- 
dures for obtaining Z s . Once this impedance matrix is evalu- 
ated, it is plugged in ( fT0| ) and the common process for solving 
the metasurface problem is carried out. The presented formu- 
lation is a general procedure allowing to solve unbiased and 
biased graphene metasurfaces under plane wave incidence. 
However, the examples considered next do not include mag- 
netic bias thus there is no non-diagonal terms due to magnetic 
field. Concerning the contribution of spatial dispersion to the 
non-diagonal conductivity terms, it is negligible in our exam- 
ples (see below) and thus the conductivity writes in essence: 



a 
cr 



(18) 



IV. APPLICATION EXAMPLES 

The remainder of the paper presents different examples of 
tunable biperiodic graphene metasurfaces, modeled and de- 
signed using the PMoM for graphene described in the previ- 
ous sections. We begin with a simple fundamental case to 
illustrate the physical basis of the considered configurations 
and their potential. Then, more realistic structures from both 
nanofabrication and biasing control point of views are consid- 
ered. Note that the effect of edges on graphene conductivity is 
also neglected, since the smallest dimensions considered are 
in the order of 1/im, which is much larger than the electron 
scattering length within the graphene layer. 



A. Freestanding metasurface 

The simplified topology, considered first, consists of a sin- 
gle biperiodic surface, whose unit cell is a cross shaped 
graphene patch (Fig. [2^ and ^p). In this example, results 
are computed assuming that the patterned graphene layer is 
free-standing in air or vacuum and the effect of gating on 
its performance is studied. The dimensions of the patch are 
L = 10 mm, d = 1.25 mm, D = 7.5 mm and a normal in- 
cidence of the plane wave on the graphene layer is assumed. 
The frequency range for the analysis is GHz < / < 30 GHz. 
Fig. [2]; shows the reflected and transmitted energy plotted in 
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terms of the excitation frequency for various bias electric field 
values. The case Eg = 20V/nm is also analyzed includ- 
ing the spatial dispersion terms, which confirms their negli- 
gible effect. The results reveal an extremely promising fea- 
ture, which is possibility to electrically control the emergence 
and strength of resonances. Note that this capability is dif- 
ferent from that presented by Ju et alP^j since there only the 
position of the resonance is controlled. Indeed, applying the 
bias electric field results in a strong increase in the surface 
conductivity, thereby exciting the resonances between adja- 
cent graphene patches. These resonances introduce maximum 
and minimum points in the reflection and transmission spec- 
tra. This behavior can be verified by visualizing the current 
profile induced on the graphene patches at / = 19.5 GHz, as 
done in Fig.|3]for three different electrostatic gatings, namely 
E = OV/nm, E = 2 V/nm and E = 20V/nm. It is ob- 
served that the induced currents gradually take on a resonating 
form when the bias electric field increases, meaning that the 
current density sharply increases at the center of the patch and 
vanishes at the edges. Though this illustration was mainly de- 
vised and presented for a first physical insight on controllable 
graphene metasurface, it is noticeable that the demonstrated 
switchable frequency selective property would find real appli- 
cations. For instance, radar absorbers can be realized, where 
the absorption peaks are switched on and off through an ex- 
ternal voltage. 

As mentioned earlier, this initial illustration is presented 
merely to introduce the advantage of graphene patterning and 
does not constitute a realistic electrically biased graphene 
metasurface. From an experimental point of view, the 
graphene layer must reside on a substrate and for simple bi- 
asing it is impractical to have DC-disconnected graphene unit 
cells. Graphene layers are sensitive to normal bias electro- 
static fields in the order of 10 9 V/m. Even relatively thin di- 
electric layer spacer can result in very large bias voltages to 
achieve the required electrostatic bias 15 . Gate voltages in the 
order of maximum 100V-200V thus require a very thin di- 
electric between graphene layer and gate, namely, in the or- 
der of some tens of nanometers. Such thin dielectrics are 
implementable and can be employed in most applications of 
graphene to realize nanoelectronic devices. However, they be- 
come problematic when considering graphene layers affecting 
the propagation of an EM wave, since a metal or silicon con- 
ductive gate layer in close proximity to the graphene layer 
would generally mask - or at least strongly affects - the de- 
sired tunable property. Therefore, the remainder of the paper 
considers metasurfaces including actual substrates, and where 
the different unit cells of a graphene layer are DC-connected 
so as to only require applying DC biasing at single location 
on the surface. The maximum DC bias voltage considered is 
200V. 

For this purpose, we propose to use periodic metasurfaces 
made of different layers of graphene 32 . Such a stack can be 
realized by first transferring a uniform graphene layer on the 
substrate and patterning it via e-beam lithography to enable 
geometrical features down to a few tens of nanometers. Then, 
the thin dielectric is deposited and finally the transferring and 
patterning steps are repeated for the second graphene layer. In 
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FIG. 2: (a) Schematic illustration of the graphene biperiodic struc- 
ture consisting of graphene cross shaped patches arranged in a two 
dimensional lattice and under a normally incident plane wave, (b) 
Unit cell of the lattice, (c) Reflected and transmitted energy versus 
frequency. The filled circles represent the results obtained when ac- 
counting for the spatial dispersion terms in the case Eo = 20 V/nm. 
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FIG. 3: Normalized induced current on the patch at uj = 19.5 GHz 
for different applied electrostatic fields are plotted for a free stand- 
ing FSS consisting of 2D lattice of cross shaped graphene patches 
exposed to a normal incident plane wave: (a) Eo = OV/nm, (b) 
E = 2 V/nm and (c) E = 20 V/nm 
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fact, the process can also be simplified by etching all graphene 
and dielectric layers in a single step. This topology has two 
main advantages: 1) It provides more freedom for tailoring the 
surface EM response and 2) it allows an efficient dynamic bias 
of the structure by applying a voltage between the graphene 
layers themselves as symbolized in Fig. [4}}. Hence, the afore- 
mentioned masking of the graphene metasurface by a neigh- 
boring metal or silicon gate is basically avoided. Obviously, 
a DC contact between the lower layer and the voltage source 
must be implemented at one location in the surface, but this 
is also very simply realized through an access cavity in the 
thin dielectric or implementing vias. Furthermore, each patch 
of the periodic metasurface must be connected to the neigh- 
boring ones for transferring the DC voltage throughout the 
whole surface. Therefore, owing to the conductive properties 
of graphene, a periodic structure with DC connected patches 
is preferred. In this case, no additional circuit or layer to sup- 
port for the DC gating is required. 



B. DC-connected double layer metasurface 



Based on the described technological approach, a graphene 
metasurface was designed to operate as a frequency selective 
surface in the far infrared range, to demonstrate the potential 
of the proposed concept at higher frequencies. The considered 
surface and unit cell are depicted in Fig.|4^ and|4j5. The same 
pattern is etched in both graphene layers, with L = 5 /xm, 
I = 0.5 fim, d = 1.25 ^m, and D = 1.5 /im. The two 
patch layers sandwich a L = 50 nm thick GaAs layer over 
an infinitely thick GaAs substrate. The relative permittivity of 
GaAs at DC is equal to e^ c — 10.8 and at far infrared fre- 
quencies measured data for the dispersive dielectric constant 
is considered 47 . Note that though GaAs might not be the best 
dielectric for graphene, in fact the performance computed here 
are weakly dependent on the dielectric permittivity. In any 
case, the results could be easily updated for other materiel 
choices, such as Parylene having similar permittivities. A 
110V bias voltage, corresponding to the normal electrostatic 
field Eq — 2.2 V/nm, is applied between the two graphene 
layers. In Fig.|4j;, the reflected energy under the excitation of 
a x-polarized plane wave is shown, when the bias voltage is 
switched on and off. The sharp resonances seen in the reflec- 
tion spectrum occurs due to the increase in the graphene con- 
ductivity. The devised graphene biperiodic layer constitute a 
controllable ultra-thin multiple band filter at infrared frequen- 
cies. The device simulation is carried out with and without 
consideration of spatial dispersion, which confirms the negli- 
gible effect of spatial dispersion in this example as well. In 
Fig. [5] the magnitude of the electric field in the xz plane is 
shown at the two different resonance frequencies. Two differ- 
ent cases of biased and unbiased graphene patches are consid- 
ered. The results evidence the considerable variation of the 
transmitted field with the electrostatic gating, which is more 
pronounced in the first resonance. 
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FIG. 4: (a) Schematic illustration of the graphene biperiodic struc- 
ture consisting of DC connected patches on both sides of a GaAs 
layer residing on an infinitely thick GaAs substrate and assumed to 
be subjected to a normally incident plane wave, (b) The unit cell of 
the lattice is depicted with the dimensions; the top figure is the top 
view of the unit cell and the bottom figure is its side-view, (c) The 
reflected energy versus frequency is plotted for unbiased case as well 
as biased electrostatic field with Eq = 2.2 V/nm. 



C. Graphene microwave absorber 

A general conclusion from the presented analyses is that 
using graphene technology ultrathin surfaces can be devised, 
which influence dramatically the incident wave. Moreover, 
in these extremely thin geometries, one has the possibility to 
electrically control the performance of the device. This further 
adds to the promise of the graphene biperiodic surfaces. Based 
on the gained insight, useful devices can be designed for spe- 
cific applications, which is the main focus of the next exam- 
ple. In Fig. [6] a multilayer absorber configuration is shown, 
which contains five layers of graphene patches separated by 
Si02 thin films with thickness t = 50 nm. The multilayer 
graphene metasurface resides on an Si02 substrate with thick- 
ness h = 3 mm backed by a perfect electric conductor (PEC). 
Fig. [(J shows the DC biasing circuit which can be used for 
applying the electrostatic bias field on the patches. 

In Fig.|6b, the magnitude of the reflection coefficient in 
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(a) (b) 

FIG. 5: Magnitude of the total electric field is shown for two res- 
onant frequencies, namely (a) 11.3THz and (b) 20.0 THz. The up- 
per figures correspond to unbiased graphene patches and the lower 
ones are obtained for an electrically biased graphene surface with 
£' = 2.2V/nm. 

terms of frequency is depicted for various bias voltages V g . 
The structure behaves as an absorber operating at /n = 
10.7 GHz when no bias voltage is applied. When a bias volt- 
age is applied, the perforated patch layers start to resonate, 
which prevents the penetration of the incident wave into the 
absorber. Thus, a large amount of the incoming field is re- 
flected, as depicted in Fig. [6};. The large variation of reflection 
coefficient at resonance frequency with the bias voltage is an 
interesting feature of the proposed absorber, which is shown 
in the inset of Fig. [6]:. According to this figure, the reflected 
energy from the absorber is controlled by the bias voltage over 
a relatively large interval. 



D. Graphene polarizer 

A remarkable advantage of graphene biperiodic structures 
is that a birefringent performance can be achieved using unit 
cells without diagonal symmetry. Based on this idea a po- 
larizer is designed as the last example, which is illustrated in 
Fig. [7] Homogeneous graphene layers are etched periodically 
to create a lattice of slots. The slots are longer along y axis 
than x axis, which leads to different performance for x and 
y polarized waves. By proper adjustment of the dimensions, 
a tunable polarizer is designed for infrared frequencies. The 
graphene patch layers shown in Fig. [7^ are fabricated on two 
sides of a Silicon Nitride (ej 30 — 6.6 and 4 R = 4.2) thin film 
with thickness t = 50 nm and the whole structure resides on a 
thick SiC>2 substrate (Fig.[7j)). The reflected energy versus fre- 
quency for V g — V and V g = 50 V are plotted for two polar- 




Frequency (THz) 
(c) 

FIG. 6: Schematic illustration of the absorber structure consisting of 
five graphene layers printed on a glass substrate: (a) The graphene 
patch layer, (b) the side-view of the absorber and the biasing circuit, 
(c) Magnitude of the reflection coefficient versus frequency for dif- 
ferent applied bias voltages (Vg) are plotted. In the subfigure, magni- 
tude of the reflection coefficient at /o = 11.2 GHz is plotted in terms 
of the applied bias voltages (V g ) for the considered absorber. 

izations in Fig.|7f;. The power reflection coefficient (R xx ) rep- 
resents the reflected energy for the electric a; -polarized wave 
(E x 7^ 0, E y = 0) when a similar plane wave illuminates 
the surface. Similarly, R yy is computed for the electric y- 
polarized reflected and incident waves (E x = 0, E y ^ 0). 
It is observed that at the resonance frequency / = 12.7 THz 
the x polarization is strongly reflected, whereas the y polar- 
ization is weakly affected. Hence, when both polarizations 
equally exist in the incident wave, the reflected plane wave 
has a dominant x polarization at this frequency. 



V. CONCLUSION 

A full-vector spectral method based on the periodic method 
of moments (PMoM) has been introduced for the analysis of 
graphene biperiodic structures. The method was then applied 
to the modeling and design of various graphene biperiodic 
structures, demonstrating that periodically patterned graphene 
structures can implement EM metasurfaces, which exhibit dif- 



8 



l|am 



50nm 



V., 




_ R with V =50V 

xx g 

_ R with V =50V 

yy a 

_ R with V =0V 

xx g 

_ R with V =0V 

yy a 



10 15 20 

Frequency (THz) 

(C) 



30 



ferent and useful dynamically-controllable capabilities, influ- 
encing both the amplitude and polarization of the EM wave. 
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